Macroscopic Quantum Tunneling of Ferromagnetic Domain Walls 



OS 
On 



O 



I 

C 

o 
o 



> 

o 
o 

o> 

-I— > 



i 

-a 
c 

o 
o 



Hans-Benjamin Braun 1 , Jordan Kyriakidis 2 , and Daniel Loss 2 
Paul Scherrer Institut, CH-5232 Villigen PSI, Switzerland 
Department of Physics and Astronomy, University of Basel, 
Klmgelbergstrasse 82, 4056 Basel, Switzerland 
(June 1997) 



Quantum tunneling of domain walls out of an impurity potential in a mesoscopic ferromagnetic 
sample is investigated. Using improved expressions for the domain wall mass and for the pinning 
potential, we find that the cross-over temperature between thermal activation and quantum tunnel- 
ing is of a different functional form than found previously. In materials like Ni or YIG, the crossover 
temperatures are around 5 mK. We also find that the WKB exponent is typically two orders of 
magnitude larger than current estimates. The sources for these discrepancies are discussed, and 
precise estimates for the transition from three-dimensional to one-dimensional magnetic behavior of 
a wire are given. The cross-over temperatures from thermal to quantum transitions and tunneling 
rates are calculated for various materials and sample sizes. 



PACS Numbers: 75.45. +j, 75.60.Ch, 75.40. Gb, 75.30.Gw 



I. INTRODUCTION 

The possibility of observing quantum mechanical be- 
havior at a mesoscopic scale has recently attracted much 
experimental and theoretical interest. First, there is the 
fundamental issue of identifying physical systems pos- 
sessing many degrees of freedom which support a col- 
lective mode that features typical quantum properties 
such as superposition behavior, interference effects, or 
tunneling through potential barriers. Well-known ex- 
amples of such systems are Josephson junctions which 
have been extensively studied in the past (for a review 
see e.g. LeggettEl). Recently the focus,-has shifted to- 
wards low-dimensional magnetic systemsa such as single- 
domain ferromagnets and antiferromagnets, but also to- 
wards nonuniform magnetic structures exhibiting domain 
(or Bloch) walls. In the latter case, one envisages a 
domain wall trapped by a magnetic pinning center — 
as provided, for example, by an impurity lowering the 
anisotropy energy locally. The domain wall can then 
escape from this potential well by tunneling through 
the energy barrier. The observability of such tunneling 
events basically depends on three conditions which can 
be stated qualitatively as follows. First, the tunneling 
barrier should be neither too high nor too wide. Second, 
the effective mass associated with the tunneling dynam- 
ics of the Bloch wall, and hence the number of spins in 
the wall, should not be too large. These two conditions 
are required in order to have a tunneling rate not too 
small, so that one can expect a tunneling event to take 
place within a reasonable amount of time (typically on 
the scale of hours or less). And third, the crossover tem- 
perature which separates the classical regime of barrier 
crossing due to thermal activations from the quantum 
regime of tunneling should realistically be in the milli- 
Kelvin range or above. Clearly, a more precise formu- 
lation of these conditions is essential since they are of 
fundamental importance for the interpretation of recent 



and future experiments in terms of macroscopic quantum 
tunneling. With this motivation, it is our goal to provide 
such quantitative estimates in the following. 

Although the idea of domain wall tunneling has, first 
been described for bulk samples in the seventies,B~L3, 
was not until the work by Stamp and collaborators! 
that this idea has received wider attention. The past few 
years have seen considerable progress in sample prepara- 
tion and has made a detailed experimental-study of the 
relaxation properties in nanowires possibleE In such ex- 
periments, the observation of a temperature independent 
relaxation (or resistance) below a critical temperature is 
taken as a strong indication for quantum tunneling. Such 
observed cross-over temperatures lie in the range of 2 
to 5 Kelvin. Resistivity measurements at low tempera- 
tures require metallic samples. The presence of conduc- 
tion electrons, however, may interfere with the tunneling 
process by providing a channel for dissipation (although 
large domain walls or low conductivities reduce this unde- 
sirable effect) . Insulating samples overcome this problem, 
but then experiments more difficult than resistance mea- 
surements are necessary. Rather than resistivity, magne- 
tization can be measured; a depinned wall which propa- 
gates down the sample will be accompanied by a sudden 
change in the magnetization. 

Theoretical estimates, based on the_same model con- 
sidered here, have been given before,0 but as we shall 
see the conclusions reached have been too optimistic. 
In particular, we find the functional dependence of 
the crossover temperature on experimentally important 
quantities such as the coercivity and domain wall mass 
are quite different from earlier calculations. This result 
has already been stated in Ref. ^, section VII, but with- 
out any details given. Below we provide these details, 
thus supporting earlier claims. The value of the crossover 
temperature is of considerable interest for the interpreta- 
tion of experimental observations since it is usually taken 
as a strong indication for the existence of quantum tun- 
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neling if the magnetization switching becomes tempera- 
ture independent.below this crossover temperature. Also, 
earlier estimatesQ predict reasonable tunneling rates for 
domain walls containing up to 10 6 spins, whereas we find 
that the number of spins in a flat wall should not exceed 
10 4 . 

The paper is organized as follows. In Sec. || we present 
the model for a ferromagnet. We then discuss the con- 
ditions under which transverse spin waves freeze out and 
the sample can be considered one-dimensional. The do- 
main wall mass is derived from the well known classical 
soliton solutions, and the origin of the impurity potential 
is discussed. In Sec. Ill we evaluate the tunneling rates 



and cross-over temperatures for a domain wall out of a 
pinning potential. Explicit numerical examples are given 
for various materials such as YIG, Ni, and in particular 
the very promising. n | er .ovskite and (badly) itinerant fer- 
romagnet SrRu03 .EKll These results are summarized in 
Tables I— III. Finally, we compare these results in Sec. IV 
with values given previously in Ref. [7| 



II. MODEL 

We consider an elongated ferromagnetic sample (or 
"nanowire") as depicted in Figure 1. We assume that 
the transverse dimensions w of the sample are small 
enough to ensure that the system behaves effectively one- 
dimensional (ID) at the typical temperatures T of an ex- 
periment. Quantitative estimates for w and T will be 
given below. Now the energy of an effectively ID ferro- 
magnet extending along the x-axis is given by 

l-L/2 

E[9, </>]=A dx\ A[{d x ef + sin 2 9(d x c/>) 2 } 

J-L/2 

+K e sm 2 e + K h sm 2 9sm 2 (j)Y (2.1) 

where A is the cross sectional area of the sample and 
the sample length L is assumed to be much larger 
than the width of a domain wall. The magnetiza- 
tion has been expressed in polar coordinates, M = 
Mo (sin 9 cos (f>, sin 8 sin 0, cos 8) with Mq = g/iBs/a 3 the 
saturation magnetization and a the lattice constant. The 
three terms in the energy density of Eq. (2.1) respectively 
describe isotropic exchange, easy-axis anisotropy (along 
z) and hard-axis anisotropy (along y). The anisotropy 
terms are of an effective nature and can contain both 
crystalline and demagnetizing (shape-induced) contribu- 
tions. A typical example of an elongated sample is shown 



in Fig|-^1 where K e 



K, 



e,cryst 



and Kh 



Kh, 



cryst 
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In the absence of dissipation, the dynamics of the mag- 
netization is governed by the familiar Landau-Lifshitz 
equations 



d t <t> = - 



7 5{E/A) 

Mn 5 COS 9 



d f cos 9 = 



7 S(E/A) 
M 8(f) 



(2.2) 



Here bot h <f> and cos 9 depend on x, t, the energy E is 
given by (|2.l|), and 7 = denotes the gyromagnetic 

ratio. 



A. ID Regime 

The system exhibits quasi ID behavior when all trans- 
verse degrees of freedom are frozen out. In order to ob- 
tain a quantitative estimate of this ID regime and thus 
of the validity of the model (^T]) , ( |2.2[ ) , we start from the 
three-dimensional (3D) spectpim of excitations around a 
Bloch wall, which is given byO 



e n ^= 2 ^JAk 2 ± +n(Ak 2 x +K e 



K h + Ak\+n{Akl + K e ), (2.3) 



and evaluate the corresponding finite size gaps. Here 
= (k y ,k z ) is the wave vector of spin waves running 
transverse to the sample, and k x is the wave vector along 
the sample. The parameter 11 = 0,1 characterizes the 
type of excitations. For n = 1 one obtains the spectrum 
of the traditional spin-waves, whereas n = leads to the 
spectrum of Winter (or flexural) wall modes which de- 
scribe a curving of the Bloch wall. In the limit of an 
infinite sample the spin waves have an anisotropy gap, 
2a 3 K e /s, while the flexural modes are gapless. For the 
finite sample widths considered here, however, the trans- 
verse spin waves and in particular the flexural modes ac- 
quire an additional finite size gap. This gap arises since 
the first excitation in transverse direction involves the fi- 
nite wave vector k mul — ir/w, where w denotes the max- 
imal transverse dimension of the sample. As a conse- 
quence of these finite size gaps, all transverse excitations 
of type n around the Bloch wall get frozen out below a 
given temperature T for sample widths w < w n , where 
w n follows from Eq. (2.3) 



W n (T) = 7T 



2A/(K h + nK e 



1 + {pTf 



nK,, 



'K h +nK e 



1/2 



(2.4) 



Here we have set p — sks/iKh + nK e )a 3 . Note that 



diverges when fc^T approaches n(2a 3 j ' s)\J K e {K e + Kh) 
(from above). 

Since the minimal energy e„ = o,k of the flexural wall 
modes is always smaller than that of the spin wave modes 
(at the same wavevectors) , quasi ID behavior at temper- 
ature T is established for sample widths w < w n= o(T). 
For instance, we find w n =o — tt-v/ A/ Kh at a "freezing" 
temperature fc^T = 2\f2Kha 3 / s, where Kh = 2ttMq for 
a slab as shown in Fig. 1. Typical numbers for Ni (see Ta- 
bic |h|) are w n=0 s» 250 A and T « IK. We note that these 
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flexural wall modes are frozen out well above the typical 
cross-over temperature T c below which one expects to see 
quantum tunneling. This cross-over temperature will be 
calculated below and is found to be of the order of 5mK. 
For a given sample width w, the transverse spin waves 
(n = 1) with k x = 0, ^ are frozen out at even higher 
temperatures. 



B. Soliton solutions and soliton mass 

In the literature, various differing values for the wall 
mass have been used. Therefore we give now a derivation 
of the wall mass from th e ex act soliton solutions!! 2 ] of 
the equations of motion ( |2.2| ). These soliton solutions 
describe Bloch walls traveling at a constant velocity v, 
and are given bycia 



%{x-vt) = 2arctane (a; - t ' t)/ l 



(2.5) 



The soliton velocity v is related to the (constant) az- 
imuthal angle 



>o by the expression, 
A jKh sin 20o 



where 



K e M y/l + K sin 2 



K h /K e 



(2.6) 



(2.7) 



We see that at finite velocities the magnetization is tilted 
out of the easy xz-plane (</>o = 0, 7r) and also that the 
Bloch wall has a limiting velocity (the "Walker limit" ) 



v [VT 



1], 



where v = 2jy / AK e /M . 

The width of the moving Bloch wall 



s = 6\i 



Ksin 2 o ] 1//2 , 



(2. 



(2.9) 



is contracted relative to the width S = \/ A / K e of a Bloch 
wall at rest. Inserting Eqs. (2.E) and (2J3) into (2.1) we 
obtain the total energy of a moving Bloch wall 



E(0 O , 



E V 1 + k sin^ 



(2.10) 



where Eq = 4Ay/ AK e is the energy of a static Bloch wall. 

If the hard-axis anisotropy energy induced by the soli- 
ton motion is much smaller than the easy-axis anisotropy, 
i.e. k sin 2 (f> <C 1, then v -C Vq*Jk, and the energy in 
Eq. (2.10) takes the form 



E{9 ,c 
with the wall mass 



M , 
E + v« 



M = A M % 



l 2 K h 



(2.11) 



(2.12) 



(provided, of course, that v < v w ). For a hard-a xis o f 
demagnetizing origin of the form Kh — 2ttMq, Eq. ( 2.12 ) 



reduces to t he Do riag expression of the wall mass Md 
(A/27:"/ 2 ) y/JQlH We emphas ize that in the presence 
of an additional strong crystalline hard-axis anisotropy, 
, with Kh cryst 3> 2ttMq, the wall 



Kt, = 2-kM, 



K, 



h, cryst ) 



mass (2.12) is substantially smaller than the Doring value 
(smaller masses lead to higher tunneling rates). Wall 
masses that are up to 10 3 smaller than the Doring value 
are fou nd in the orthoferrites.Ej 

Eq. (2.11) shows that a moving domain wall behaves 
like a particle of mass M. The dynamics of the domain 
wall 6q(x — X) with X = X(t) can therefore be described 
by the action of a free particle of mass M, 



5 (0) 



M ■ , 

at -x 2 . 



(2.13) 



For a microscopic derivation of Eq. ( 2.12 ) from the quan- 
tum spin action within a coherent spin state path integral 
formalism and a collective coordinate technique (and also 
including the effects of dissipation via spin waves) we re- 
fer the reader to Refs. [lqjF7|,H- 



C. Impurities and pinning potentials 

So far we have focused on an ideal sample with per- 
fect translational invariancc. In realistic samples this in- 
variance is broken by impurities or modulations of the 
sample cross section. We extend now the above consid- 
erations to this situation and discuss the effects of an 
external magnetic field. For simplicity we treat first a 
point-like impurity, consisting of a single atom at x = 
with easy-axis anisotropy K p ^ K e . Such an impurity 
can be described by changing the anisotropy constant in 
Eq. (|2.lD in the following way, 



K e ^ K e + K p {*), K p (x) = -V&tf(x). (2.14) 

where Vq = (K e — X p )a 3 , for v such impurities we ev- 
idently have Vo = v(K e — K p )a 3 . Without loss of gen- 
erality, we consider in the sequel the case of attractive 
impurities, i.e. Vq > 0. 

A uniform external field along the z (easy)-axis is de- 
scrib ed b y a Zeeman term —HMq cos 9 in the energy den- 
sity © 



Both pinning and external field thus lead to 



the additional energy 



E' = J d 3 x^K p {x)sm 2 9 - HM a cos6^. (2.15) 

The impurity now breaks the translational symmetry 
transverse to the sample. We consider here the situa- 
tion of weak pinning where the pinning energy is much 
smaller than the static wall energy, Vq/Eq <S 1. In this 
case deviations from the flat Bloch wall configuration 
(#o, (j>o) can be neglected. Note that Eq = 2NK e a 3 where 
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N = 2A8/a? is the number of spins in the static wall. 
Therefore the weak pinning assumption can be satisfied 
even in the case of many impurities as long as the con- 
centration of impurities within the wall volume is small. 

To lowest order in Vq/Eq, we can then insert the static 
soliton solution QJrf — X) into Eq. ( 2.15 ) and obtain the 
additional energyEil 



E 1 = V P (X) - hX 
with the pinning potential 

V P {X) = -Vq sech 2 



X 



and the force due to the external field 
h = 2AM H. 



(2.16) 



(2.17) 



(2.18) 



In ( |2.17| ) we have used that 6 = S(1 + O(v/v ^/k) 2 ). Note 
that even a point-like pinning center of the form (2.14) 
creates a shallow potential ( p. 17 ) varying on the length 
scale of the Bloch wall wi dth 5 = -J ' A/K e . 

The pinning potential (2.17) not only holds for point- 
like impurities but also describes pinning due to varia- 
tions in the cross-section if they extend over length scales 
I shorter than the domain wall width 5. To be specific, 
let us consider a constriction where the cross sectional 
area A(x) = A — AA(x) is locally reduced, i.e. AA(x) 
vanishes for \x\ > I. Let Av — JdxAA(x) denote the 
missing sample volume of the constriction. The total 
wall energy is then 



E' = Jdx A{x) {2K e sin 2 6 - HM cos 9} 
2 X 

= —2K„ Av sech — — hX + const. 



(2.19) 



where 8 — 8q(x — X). Thus the e ffe ct of the constriction 
is again described by Eqs Q2.16Q , fl2.17| ) but now with 
Vq = 2K e Av = EqAv/2A5. The weak-pinning limit is 
thus justified as long as the volume Av is small compared 
to the volume 2A5 occupied by the domain wall. In the 
second line of Eq. (2. IE), we have suppressed a small 



additional Zeeman term (hAv/2A) tanh(Jf/<5) which is 
an irrelevant constant for large X, while for small X it 
renormalizes ft by a factor 1 — Av/2A5. However, this 
renormalization is small in the weak pinning limit con- 
sidered here and thus experimentally not relevant. 

In conclusion, we find that the dynamics of a domain 
wall in an external field, and in the presence of a point- 
like impurity (or a constricted cross section) is described 
by the action S = — J dtE' . Explicitly we have 



S 




V P (X) + hX 



(2.20) 



with V p as in ( [2.17 ). Vq depends on the impurity or 
constriction parameters as defined above. The dynamics 
of the Bloch wall is thus seen to be equivalent to that of 
a particle of mass M in a potential V p under a force h. 



III. DEPINNING VIA QUANTUM TUNNELING 

In this section we calculate the tunneling probability 
of a Bloch wall out of a pinning potential V(X). For the 
moment, let us consider a pinning potential of arbitrary 
shape, as might arise e.g. from the presence of many ran- 
domly located impurities. We shall return below to the 
specific case of the generic sech 2 potential. 

Interested in tunneling phen omen a, we consider the 
Euclidean version of the action ( 2.20| ) 



S E 




f) — > 



where units have been chosen such that h 
potential energy for the wall is given by 

U{X) = V(X) - hX. 



(3.1) 



1. The 



(3.2) 



In E q. ( | 3.1[) , /3 = the wall mass M is given by 

Eq. Q2.12j ), and V(X) is some smooth pinning potential 
which, for the present, we keep arbitrary. It is only 
assumed that V(X) < and that it tends to zero for 
\X\ — > oo. It then follows that V(X) has at least two 
inflection points Xi, defined by V"(Xi) = 0. 

Let us consider the situation of a vanishing external 
field where the wall is pinned at a local minimum X m j n 
of V(X). Let Xi be the inflection point closest to the 
right of X min . Thus V'(Xi) > and V&(Xi) < 0. 
At small (positive) values of the external field, the wall 
is still trapped at X m j n , but as the field is increased, 
the potential becomes increasingly tilted and finally, the 
metastable state ceases to exist at the coercive force 



(3.3) 



where h c = 2AMqH c > 0, with H c the classical (zero 
temperature) coercivity. In Fig. ^| we plot the poten- 
tial energ y U(X) with the sech 2 pinning potential of 
Eq. (2.17). The three curves shown are for external fields 
near the classical coercivity. It should be kept in mind, 
however, that the following analysis is valid for arbitrary 
pinning potentials (subject to the conditions expounded 
in the preceding paragraph). 

The possibility of quantum tunneling arises when the 
external field is close to the classical coercivity, i.e. 

< e = 1 - H/H c < 1. (3.4) 

The potential ( |3.2| ) can then be expanded around the 
inflection point Xi of U (X) to yield 

U(x) = if< s >(X,)x» + (-i^V^X,))" 2 * 2 

Several comments are in order regarding these expres- 
sions. First, the third derivative V^(Xi) < in general 
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depends on the coercivity. Also, we have shifted coor- 
dinates so that the minimum is now U(0) = 0. In the 
second line, we have introduced the tunneling distance 
d > 0, defined by U(d) = 0, aad the barrier height V, nax . 
These are explicitly given byt3 



2 5/2 (/j c e)3/2 



-VW(Xi)]V*' 



3V2- 



(h c e) 



1/2 



-V^iXi)} 1 / 2 ' 



(3.6) 



(3.7) 



For external fields close to the coercivity, the Euclidean 
action associated with the tunneling of the domain wall 
is thus given by 



S E [x] = [ 
Jo 



"Jut* 

2 Ut 



U(x) 



(3.8) 



with U(x) as in Eq. (3.5). This action is rendered sta- 
tionary by the "bounce" trajectory 



x b ( T ) = d sech 2 ujbT 



(3.9) 



which runs from x = to x = d and back to x = 
for t increasing from —0/2 to (3/2 with (3 — * oo. The 
characteristic tunnel frequency is given by 



LO b = (3/2) 3 /VKnax/A« 2 



(3.10) 



Note that u>b is half the harmonic oscillation frequency 
in the potential minimum of U. The tunneling action 
So = S[xb] can be calculated without explicit knowledge 
of the above bounce trajectory, i.e. 



For the observation of quantum tunneling it is also 
important to ensure that the thermally activated tran- 
sition rate over the barrier, Pr — u>o ex p[ — V ma ,x/kBT], 
does not exceed the tunneling rate P through the barrier. 
This is the case if the sample temperature T is less than 
the cross-over temperature T c , which can be estimated 
by equating the corresponding transition rates. By as- 
suming that the prefactors are approximately equal we 
have (after reinstating h) the relation ksT c = V max h/So, 
which yields 



knT r 



5 

8d 



2V m 



3M 



18 



-TlLOh 



(3.14) 



In order to obtain further quantitative understanding, 
we now apply the above results to the specific case of the 
generic pinning potential V p (X) = -V sech 2 f . With 
Eqs. (3.3), (|3.6|), and (3.7) we immediately find that 



4 Vp 
3 3 / 2 S ' 



V m 



2V2 



36 



(3.15) 



The coercive force is thus linked to intrinsic properties 
of the pinning potential — the ratio of potential strength 
Vq and characteristic length scal e 5. Comparison of these 
expressions with Eqs. ( |3.7| ) and (3.6) shows now explicitly 
that V^(Xi) = —4h c /r indeed depends on the coercive 
force. 

The tunneling exponent, cross-over temperature, and 



tunneling frequency follow from Eqs. ( 3.1 2| ) , (3.14), and 
( |3.10| ) and are given by 



S = (6/5)hNsy/H c /2nMo (2e) 5/4 , 



(3.16) 




So = 2V2M / dx y/U{x), 



dy/MV max = 



18 V ma) 
5 wb 



(3.11) 
(3.12) 



Note that the factor of 2 in the first equation arises be- 
cause the escape rate is determined by the action over 
the whole bounce which leads from x — to x — d and 
back to x = 0. 

T he e scape tunneling rate P for the potential U in 
Eq. (3.2) has been, calculated before in a different context 



[see e.g., Weissjia p. 109, Eqs. (8.12) and (8.16)]. It is 
explicitly given by the standard WKB expression 



P = 4cj b V 1 55 /27re"' So . 



(3.13) 



Typically, quantum tunneling will be observable if the 
time between tunneling events, i.e. the inverse escape 
rate -P -1 , does not exceed a few hours. For a typical 
attempt frequency cob (given approximately by the expo- 
nential prefactor in P) of the order of 10 9 s _1 this requires 
that the exponent Sq/% be less than about 30. 



k B T c = (5/18) ff /W27rM tf c (2e) 



1/4 



uj b = (g^ B /ny2^MoH c {2e) 1 '\ 



(3.17) 



(3.18) 



where N = 2A8/a? is the total number of spins in the 
wall and we have assumed a purely shape-induced hard- 
axis anisotropy, i.e. — 2ttMq for a slab geometry. 
Alternative but equivalent expressions for the bounce fre- 
quency u>b, the WKB exponent Sojh, and the crossover 
temperature T c are listed in Table m. 

To illustrate the above analytical results with concrete 
numbers we have collected in Tables |n| and |T3 various 



values for several ferromagnetic samples of the shape 
shown in Fig. 1, namely Yttditm Iron Garnet (YIG), 
Nickel, the perovskite SrRu03,c2rLil and "large easy-axis" 
materials considered in Ref. 0. From Table [II it becomes 
clear that the typical number N of spins one can expect 
to tunnel coherently out of a pinning potential within 
reasonable time (a few seconds) is of the order of 10 4 or 
less, and that the associated cross-over temperature T c 
is typically less than lOmK. A stark exception to this 
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is SrRu03, in which our theory predicts 10 3 spins can 
coherently tunnel approximately once every millisecond 
with a crossover temperature of 37 milli-Kelvin. It would 
therefore be very interesting to look for domain wall tun- 
neling in this mate ria l. 

From Eqs. (|3.16|) -( ^18|) we see that in order to opti- 
mize the observability of quantum tunneling it would be 
desirable to have materials that possess both a large co- 
ercivity H c and a large hard axis anisotropy Kh but with 
the ratio Hc/^/K^ being small. Such materials would 
have a small WKB exponent (i.e. a high tunneling rate) 
and a high cross-over temperature. There is, however, 
some leeway by carefully choosing the experimentally 
tunable parameters N and e (see e.g. SrRu03 in Tables |l| 
and HI . 



IV. DISCUSSION 
A. Comparison with earlier results 

A discussion of domain walls tunneling out of impurity 
potentials has been given by Chudnovsky, Iglesias and 
StampQ which we shall henceforth refer to as CIS (for 
earlier work in this field see also the references contained 
in CIS). Our results presented in Table III for experi- 
mentally vital quantities such as the WKB exponent Sq 
and the cross-over temperature T c to the quantum regime 
differ substantially from the results given in section VII 
of CIS. Therefore a comparison of the two approaches 
appears necessary. 

Before turning to the most crucial difference between 
the two approaches — the functional dependence of the 
pinning potential on the coercivity — let us first remark 
that the mass M CIS used in CIS differs from the soliton 



mass M given in Eq. ( 2.12 ). In fact 



M cis = M _ 



(4.1) 



The two masses agree only in the limit of large hard- 
axis anisotropy, i.e. k = Kh/K e 3> 1. However, in 
the experimentally important limit of k <C 1 we have 
M CIS /M = 4/k ^> 1. Thus, we would expect our action 
to be smaller and the cross-over temperature larger than 
the CIS results. 

However, this tendency is more than compensated by 
the CIS assumption that the pinning potential width w 
is independent of the coercivity. For a field close to the 
coercivity, CIS find (their Eq. (25)) for the total pinning 
potential of arbitrary shape 



A W U(Z) = ctqA, 



fV2h 



CIS, 



2w 



— (- 

6w 2 \ S 



(4.2) 



where = Sh c /Eo, (JqA w = Eq, and A w — A. w is 
a parameter which is assumed to be independent of h c 



and is set equal to unity i n se ction VII of CIS. However, 
comparison with our Eq. ( |3.5|) reveals that the parameter 
w 2 = —aoAw/V^ (Xi)S 3 is not arbitrary but depends on 
the details of the pinning potential. 

In particular, for the generic sech 2 pinning potential 
( |2.17 ) it follows that w is coercivity-dependent: 



CIS 



(4.3) 



Since experimentally h~r = H c /H a (where H a = 
2K e /Mo) lies typically in the range 10 _2 -10 -5 , this im- 
plies a value of w in the range 5-100. 

By assuming w ~ 1 in their final section, CIS have 
thus a priori fixed the depth of the pinning potential 
to the extremely large value of a third of the total wall 
energy, Vo = Eq/3. This situation corresponds to a re- 
gion of length 4(5/3 extending across the entire cross sec- 
tional area being replaced with magnetic ions of vanish- 
ing anisotropy. 

If we insert our expressions for M and w into Eq. (27) 
of CIS, then their tunneling action agrees with ours (up 
to a minor error of a missing factor of 2 1 / 4 in Eq. (27) 
of CIS). With these substitutions we also find agreement 
between our T c and Eq. (44) of CIS (where again a factor 
2 -1 / 4 is missing). 

Our results strikingly differ from CIS when it comes 
to the explicit computation of experimentally important 
quantities. In Eqs. (88) and (87) of CIS, numerical fac- 
tors 2 1 / 4 48/5 ~ 11 and 5\/2/2 1 / 4 36 ~ 1/6 respectively 
are suppressed compared to CIS Eqs. (27) and (44) for 
iSo and T c . Together with their assumption that w ~ 1, 
this leads (for a material with h^ ls = 10~ 3 ) to an under- 
estimation of the WKB exponent Sq by a factor of 700 
for a planar domain wall tunneling through the poten- 
tial of a single defect. At the same time the cross-over 
temperature T c between quantum tunneling and thermal 
depinning is overestimated by a factor of 24. 

Finally, we mention that while our cross-over temper- 
atures and WKB exponents differ substantially from the 
values presented in CIS, we find the same scaling of T c 
and Sq with respect to the reduced field e. 



B. A mechanism for increasing T c 

Explicit expressions for various materials such as Ni, 
Fe, YIG, and SrRu03 are presented in Table [II. It 
is seen that the transition temperatures are in the mK 
range. A notable exception is SrRu03, whose small do- 
main wall width leads to a narrow potential well and thus 
to a considerably higher transition temperature — around 
40mK. 

We must mention, however, that our simple analysis 
here does not rule out discernible tunneling of larger walls 
at higher temperatures. Recent experiments on domain 



G 



wall dynamicsu have been interpreted as evidence of tun- 
neling. This evidence primarily stems from the occur- 
rence of temperature independent phenomena below a 
cross-over temperature of a few Kelvin — three orders of 
magnitude larger than our estimates here. Rather than 
comment on these experiments, we will instead discuss a 
plausible mechanism which may raise T c and/or decrease 
So. 



So = S ok V(l + m)/A 3 



(4.9) 



Table III shows that the tunneling rate becomes ap- 
preciable if the tunneling distance is smaller than 10A. 
Thus, variations of the pinning potential V(X) on this 
length scale could dramatically affect the tunneling be- 
havior. Such variations do not occur for a random su- 
perposition of pinning potentials of the sech 2 type. How- 
ever, under certain circumstances the underlying crystal 
lattice can provide such a modulation, in particular if the 
wall width is only a few lattice constants. 

To get some quantitative ideas about the consequences 
of a modulation with the period of the lattice constant, 
let us add the term V VCI (X) = V\ sin(27rA/a + Q) to the 
pinning potential in Eq. ( fL2| ) Here, £ is a phase which 
we conveniently choose as C, — ~2irXi/a, where AQ is the 
inflection point of the sech 2 potential, i.e. sech 2 Xi/5 = 
2/3. Thus Xi remains the inflection point of V + V pcl . 
This new potential has associated with it a coercivity 
whic h is a factor of (1 + /i) larger than the coercivity in 
Eq. (|3T5l), i.e. h c = hf d (l + fi), where 



M = 



27 7T ViS 
Vo~a 



(4.4) 



In the vicinity of Xi, we can carry out the same expan- 
sion as outlined above to again obtain the cubic potential 
given by 



V(x) = V2e^r\x 2 - ^AV. 
6 3o z 



(4.5) 



The new parameter A, which equals unity for V\ = 0, is 
given by 



1 + r 2 fi 
1 + A* 



1/2 



(4.6) 



where r = nS/a. In Eq. ( |4. 5| ) we have also redefined the 
reduced field e to reflect the increased coercivity, i.e. 



! h _ e old + /i 
h c 1 + /i 



(4.7) 



Note that we now remain in the tunneling regime (i.e. e > 
0) even if e old becomes negative. Eq. (4.5) gives rise to 



a new tunneling distance d and also a new barrier height 
Umax, which results in a larger crossover temperature T c 
and a reduced WKB exponent So ■ These new expressions 
are explicitly given by 



wher e T ° and Sq contain the new definition of e, 
Eq. (4.7), since it is now e (and not e old ) which is the 



T c = T C °V(1+M)A 



(4.8) 



experimentally tunable parameter controlling the barrier 
height for tunneling. To obtain a numerical estimate of 
this effect, we need to estimate the magnitude V\ of the 
periodic piece relative to the impurity strength Vq. Let 
us take, for example, V\ = 10 _1 Vo and a domain wall of 
ten lattice constants, S = 10a. In this case, T c « 10T° ld 
and So « l(T 2 S£ ld . The estimates here for the WKB 
exponent must be taken with caution because these are 
estimates for the tunneling through only one of the (pe- 
riodic) barriers. Due to the shape of the impurity poten- 
tial, one should expect that there are ~ 8 / a such barriers 
to tunnel through before the wall is free. Assuming inco- 
herent sequential tunneling, So will effectively increase 
by a factor of about 10 for the estimates just given, 
i.e. Sq S « lOiSo- Nevertheless, this very simple argu- 
ment shows that observation of tunneling of larger walls 
at higher temperatures is not necessarily ruled out. 

The case of a periodic potential superimposed on the 
impurity potential is interesting from another perspec- 
tive. For e old = (but e > 0), the soliton originally 
pinned at the impurity now sees an effectively flat peri- 
odic potential through which it may tunnel coherently. 
Since the soliton is a particle- like excitation, we basically 
have here the physics of a particle in a periodic poten- 
tial, and thus the seliton can form Bloch-like states of 
definite momentum!] Increasing the field still further in- 
duces a force on the soliton and the possibility arises of 
Bloch oscillations of a magnetic soliton. This idea was 
first laid out semiclassically in Ref. [l6] and essentially the 
same physics, holds down to the extreme quantum case 
of spin- 1/2 E3 

Finally, we briefly mention that thermally assisted tun- 
neling may also raise the effective T c . If the pinning po- 
tential contains some internal level structure, then the 
wall may be thermally excited to some higher level, and 
only then tunnel out of the pinning site. A more detailed 
analysis of this problem (e.g. along the lines discussed 
in the context of MQT in SQUIDS, see WeissN) is re- 
quired, however, before concrete statements concerning 
T r can be made. 



V. CONCLUSIONS 

We have given in this paper a detailed derivation of the 
tunneling problem of a planar Bloch wall out of a pinning 
potential. We have focused exclusively on a quasi-one- 
dimensional ferromagnet with biaxial anisotropy, and 
have given estimates on when a system can be consid- 
ered quasi-one-dimensional. In particular, the flexural 
spin-wave modes, while gapless in infinite systems, ac- 
quire a gap for the finite geometries shown in Figure 1. 
If the sample temperature is below this energy gap, then 
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the flexural modes cannot be excited and can hence be ne- 
glected. For cross-sectional areas on the order of 10 4 A 2 , 
this energy gap is much larger than the crossover tem- 
perature at which quantum and thermal transitions are 
equal. 

We have modeled the pinning center as an impurity 
which decreases the easy-axis anisotropy at a single point 
in space. Even such a point-like impurity produces a shal- 
low pinning potential which varies on the length scale of 
the Bloch wall width. We have related both the height 
and width of this potential to a coercivity. Detectable 
tunneling can only occur if the external field is very close 
to this coercivity, i.e. we must have e = 1 — H/H c on 
the order of 10~ 2 or 1CP 3 . This is an important number 
to determine experimentally. For example, if an experi- 
ment has e < 0, then one would also observe temperature 
independent depinning, but of course this cannot be as- 
cribed to quantum tunneling — it is trivially due to the 
vanishing of the barrier height. 

Within an instanton approach, the WKB exponent, 
bounce frequency, tunneling rate, and crossover temper- 
ature have been calculated and different analytical forms 
for these quantities can be found in Table GL We have also 
given estimates of these quantities for specific materials. 
The material parameters can be found in Table |J, an d 
the estimates in Table [II. In particular, the perovskite 
SrRuOa seems a promising candidate and we hope this 
work can motivate some experimental studies into this 
material. _ 

We have compared our results with previous world] and 
have found our calculations to predict lower crossover 
temperatures and a lower maximum number of spins 
which can coherently tunnel out of a pinning potential. 

Finally, we have briefly discussed how the effects of a 
periodic potential, may lead to a larger crossover tem- 
perature and a smaller WKB exponent. 
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x=X 



FIG. 1. Shown is a Bloch wall configuration in a thin long 
slab, i.e. 8q(x — X), and <j>o — n. 
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FIG. 2. We plot the potential energ y U( X) of Eq. ( |3.2| ) 
for the sech 2 pinning potential of Eq. (2.17). For this pin- 



ning potential, the coercivity h c equals (4/v27)Vb/<5. The 
three curves show the potential energy for the external field h 
slightly above, right at, and slightly below the coercivity field 
h c . By expanding about the inflection point (shown by small 
circles in the figure) , these curves are very well approximated 
by a cubic potential, as discussed in the text. 
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TABLE I. Summary of equivalent expressions for bounce frequency u>t,, WKB exponent So/ft, and cross-over temperature 
T c . N — 2SA/a s denotes the number of spins in the wall. In the last column it has been assumed that Ku = 2nM§, i.e. the 
sample has the slab geometry of Fig. 1 and there is no crystalline hard-axis anisotropy, ifh.cryst = 0. 





f3\ 3 / 2 1 /V milx 
V 2 J d V M 


V MS \ 2 / 


(f) V4 


2 3 / 4 9M^M 05 re 1/4 
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2 3/4 f^V^ £5/4 


T c 


/ 2 5 h 1 / V' m „ 
V 3 8 k B d V M 


o3/4 5 h /~~hT 1/4 
36 fc B V MJ e 


18 k B V Af W 





TABLE II. Saturation magnetization Mo, shape anisotropy = 2-kMq for a thin film, easy-axis anisotropy constant K e , 
exchange A, wall width S = yj A/ K e , wall mass M, and coercivity _H" C for various materials. 





Mo 
[Oe] 


[10 5 ^1 


K e 
[10 5 -^] 
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[A] 


M/area 
[lCn 10 -IU 

rm - -1 


H c 
[Oe] 


YIG 
Ni 

large 

K e e 


196 a 
508 c 
200 


2.4 
16 
2.5 


0.25 b 
8 d 
100 


0.43 a 
1 
1 


414 
112 
32 


1.2 
4.6 
16 


10 
100 
10 


SrRuOs 1 


159 


1.6 


20 


0.023 


11 


48 


10 4 



a Ref. 23, p. 65 
b Ref. 24, p. 313 
c Ref. 25 p. 270 
d Ref. |f p. 569 

e Example given in Ref. |^, parameters taken from there except for A which has been replaced by the most common value 
f Ref. |g 



TABLE III. Cross sectional area A, number of spins TV in the wall, e, tunneling distance d, cross-over temperature T c , WKB 
exponent So/Ti, oscillation frequency ui, and inverse tunneling rate P _1 for various materials. 





A 


e 


d 


T 


So/h 




p-i 




[A 2 ] 




[A] 


[mK] 




[io 9 ^- 1 ] 


M 


YIG 


50 x 200 


io- 1 


280 


3 


1268 


2.6 


CO 




[TV = 3.4 ■ 10 4 ] 


10 -2 


88 


1.6 


71 


1.5 


2 ■ 10 20 






5.7- 10~ 3 


66 


1.4 


31.1 


1.3 


1433 






10~ 4 


8.8 


0.5 


0.2 


0.47 


2 • 10 -9 


Ni 


50 x 200 


io- 1 


75 


14 


1740 


13 


oo 




[TV = 2.4 ■ 10 4 ] 


10~ 2 


23 


8 


98 


7.6 


10 31 






3.9- 10~ 3 


15 


6.3 


31.1 


6 


310 






io- 4 


2.4 


2.5 


0.3 


2.4 


3-10" 10 


large K e 


50 x 200 


io- 1 


21 


3 


98 


2.6 


5 • 10 31 




[TV = 8.0 x 10 4 ] 


3.6- 10" 2 


13 


2.1 


31.1 


2 


931 






10~ 2 


7 


1.6 


5.5 


1.5 


2 • 10~ 8 






10~ 4 


0.7 


0.5 


2-10~ 2 


0.47 


5 • IO" 9 


SrRuOs 


50 x 200 


io- 1 


7.4 


79 


873 


37 


CO 




[TV = 3.4 x 10 3 ] 


10~ 2 


2.3 


44 


49 


21 


10 9 






5 x 10~ 3 


1.7 


37 


21 


18 


10" 3 






IO -3 


0.74 


25 


2.8 


12 


10 -io 
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